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SUMMARY 


A plane elastic prohlem of bodies bounded by concentric circles and heated 
at the inner boundary is solved in the "uncoupled, quasi-static" approximation. 
The heating takes place by surface heat transfer, with limiting cases corre- 
sponding to sudden temperature change and to constant flux of heat. The stress 
boundary conditions treated are those of constant pressure at the inner wall 
while the outer wall may be either free or rigidly restrained. 

The time -dependent results for stresses, displacement, and temperature are 
calculated and shown for several cases with varying values for the parameters 
specifying heat transfer and wall thickness. The effect of varying Poisson* s 
ratio is also examined. Formulas for steady- state distributions of the quanti- 
ties of interest are also given, as well as some approximate results that can 
be used for very small times. 


INTRODUCTION 


There are many engineering uses for cylindrical vessels or annular disks, 
and some of these involve differential heating of inner and outer cylindrical 
surfaces with consequent thermal stresses. Some examples are given by Heisler 
in reference 1. The present work was begun in an effort to calculate the 
thermoelastic behavior of some disks used in a constricted arc wind tunnel, in 
which arc-heated gas flows through the central hole of the disks. The formula- 
tion of the problem and the subsequent solutions are, however, cast in general 
terms and the results may be adapted to a variety of specific applications. 

The calculation of thermoelastic stresses in circular cylinders has 

received a good deal of attention in the past, as attested by the contents 

and/or bibliographies of references 1-8. Most of the earlier work was con- 
cerned with the solid cylinder (or disk) and the paper of Jaeger (ref. 9) is a 
useful collection of numerical and graphical results for that case. Hollow cyl- 
inders (or annular disks) have not, however, been so well covered. The early 
papers of Dahl (ref. 4) and Kent (ref. 7) give analytical results for the pres- 
ent problem, but in reference 4 the numerical results apply to either thin- 
walled cylinders or to thick-walled ones with special temperature distributions. 
In reference 7j there are some very brief numerical examples given. 

More recently, Trostel (ref. 10) has given a general solution for a hollow 
circular cylinder of finite length, from which are extracted some results for 

the infinite-length, or two-dimensional, case. A method for approximate solu- 

tion of the (two-dimensional) hollow cylinder problem (among others) is set out 
in reference 11, and the case of a thin-walled cylinder with internal heat gen- 
eration is treated in reference 12. This list is certainly not exhaustive, but 
serves to illustrate the fact that there does not seem to be any systematic 



presentation of thermoelastic stresses for the two-dimensional body bounded by 
concentric circles and undergoing various types of heating. This is probably 
because, as noted in reference 1, there is a new parameter to be considered 
over those for the solid body, and the resulting analytical expressions are 
rather formidable unless there is an electronic computer available. 

It was therefore considered that an extensive set of calculations of ther- 
mal stresses in the two-dimensional elastic case (annular disks or long hollow 
cylinders) would be useful in several fields of engineering. This is particu- 
larly so for the disks with thick walls where approximate results are not so 
simply derived as for thin walls. The theory to be used in solving this prob- 
lem is the one aptly called "quasi-static, uncoupled" in reference 2. In this, 
the dynamic effects of heating, and the effect of volume change on heating, are 
ignored. As a result of this approximation, the heat-conduction problem can be 
solved first, and the result used in determining stress. Circumstances under 
which this decoupling is permissible are discussed in reference 2. There, 
examples in which more exact solutions are obtainable are analyzed, and some 
indications as to the size of the error are deduced. 


In a rough way, it appears that the quasi- static uncoupled solutions are 
valid when the product of the combinations 
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(Tj_ is the initial, strain-free temperature, p the density, c v the specific 
heat at constant volume of the material, e the rate of volume dilatation, and 
T the rate of change of temperature. The remaining symbols are defined in the 
appendix.) The smallness of the second ratio means that pronounced lag between 
temperature and displacement should not arise, nor should vibrations be induced. 
(See ref. 2 for a complete discussion.) 


The manner of heating the body in question is of some importance. In much 
of the thermal stress work the initial condition of sudden change of tempera- 
ture at some surface is taken. The overly conservative design requirements 
from this condition are pointed out by Nadai (ref. 13 , pp • 397 “ 398 ) and by 
Heisler (ref. l) . It is more realistic to consider that there is a surface 
heat-transfer resistance, and to alter the boundary conditions accordingly. 

This causes some increase in analytic complexity and, as a result, a new param- 
eter appears, namely the Nusselt, or Biot, number which involves the surface 
heat-transfer coefficient. 


Once the transient temperature distribution is known, the stresses may be 
determined. For the case studied, this amounts to evaluating certain integrals 
of the temperature distribution (as given, e.g., in ref. 1^) together with con- 
stants to fit stress boundary conditions. These conditions are here taken to 
be that there is a fixed pressure acting on the inner surface, while the outer 
surface may be either free or rigidly restrained. Intermediate cases of elas- 
tic restraint can be interpolated between these extremes. 

Formulas will be given from which temperature, displacement, and stresses 
can be calculated for arbitrary position and time. Numerical results which 
sample the full range of the parameters affecting the character of the solution 
will be presented and discussed. These should be sufficient for a designer to 
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obtain estimates for stresses in cases other than those directly shown* Steady- 
state results are quite simple and are easily calculated for any case of interest* 

At the opposite extreme from the steady- state distributions are the ini- 
tial ones, that is, the temperature and stress distributions which occur 
shortly after the onset of heating (or cooling) • The starting distributions 
are of considerable interest because the gradients of temperature and hence of 
stress are most severe then (when there is a surface heat-transfer resistance)* 
They are difficult to calculate by use of the general results, owing to slow 
convergence of series expansions. However, operational techniques (see 
Goldstein, ref. 15 ) lend themselves to approximate calculation of these short- 
time results, and a number of such formulas are given. 

ANALYSIS 


As stated in the Introduction, the quasi- static , uncoupled approximation 
in thermoelasticity will be employed* In this, the temperature distribution is 
first determined as a solution to the linear heat- conduction equation, and then 
the temperature function is used to determine displacement and stresses. 

Temperature Distributions 

First, the temperature distribution for a region bounded by concentric 
circular cylindrical surfaces without variation along the axial direction can 
be written out* The problem is also specialized to allow only radial variation 
of all quantities. There is to be a choice of three boundary conditions on the 
temperature or its gradient (flux) at the inner cylindrical surface. The equa- 
tion of heat conduction is, in the absence of extraneous heat sources, 

5 2 t + l = i St 

Sr 2 r Sr K St 

where the symbols are defined in the appendix. The thermal diffusivity k is 
taken to be constant. It will be convenient to introduce the dimensionless 
variables (see sketch) 

P = a" 1 * 1 , P 0 = a' 1 !: , t = Ka _2 t 

where a is the outer radius of the 
disk and b the inner. The equation 
of c onduc t i on bee ome s 


SfT + i _ St 

Bp 2 P C*p 5 X 

As initial condition, always 

T(p,0) = T x 


( 1 ) 


( 2 ) 
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•where T x is a temperature at which the material is strain free. A boundary 
condition at the outer surface r = a, hence at p = 1, is always to he 


T(1,t) = T, 


( 3 ) 


The boundary condition at the inner surface, r = b, or p = p Q , will be one 
of three : 


Case A 


t (p 0 j t ) -^o 


Case B 

-k|S 

= aF Q 


dp 

P=Pq ° 

where F 

o 

is the heat flux. 

Case C 

-k — 

= ha| T 0 ■ 


dp 

P=Po L 

or 




St 

m frr 


dp 

P=Pq P °l- 


Tq “ T( p Q , t) 


(4a) 

(4b) 


(4c) 


where h is the surface heat-transfer coefficient and m (= hb/k) is the Biot 
number. The boundary conditions are here taken to be independent of time, for 
t > 0. More complicated time variations can be treated by use of Duhamel ! s 
theorem (ref.16) or by direct use of the convolution process for Laplace trans- 
forms. In the last boundary condition, C, let m -» oo . Then the flux may be 
loosely considered as becoming infinite such that the temperature at p = p Q 
jumps suddenly to T Q at t = 0, and thereafter remains at T 0 # Thus, this limit 
corresponds to case A. On the other hand, if the boundary conductance h -* 0 
while the reservoir temperature T q becomes very large, it can be thought, 
again rather loosely, that this corresponds to case B, with constant flux. 

Thus, when a solution for case C has been obtained, the two limiting cases 
m = oo, o should lead to solutions for cases A and B, respectively. This can, 
of course, be checked by actual calculation. 

The problem posed in equations (l) to (4) can be solved by standard tech- 
niques; in fact, a solution for a more general problem, from which the present 
one can be derived, is given in reference lo. However, in view of the need 
which will arise later for the Laplace transform of the solution, it is conve- 
nient to go very briefly through the process. The Laplace transform of a func- 
tion is here denoted by a bar, thus 

pOO 

f(p) = / e -PT f (OdT 

oo 
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In terms of these new variables T and p ; it is now necessary to solve the 
differential equation 


where 


afr 

dp 2 


+ i — - q 2 T = -T, 
P dp H 1 


(la) 


q 2 = p 


and to make the solution conform to the appropriate boundary conditions at 
p = p 0 ,l. The last equation is recognized as Bessel's for modified functions 
of order zero. The complete solution is then 


T = AI 0 (qp) + BK 0 (qp) + p _1 T 1 


where I 0 , K 0 denote the modified Bessel functions of first and second kind as 
defined by Watson (ref. 17)* The determination of the constants A, B is made 
by means of boundary conditions; only condition C (eq. (4c)) will be consid- 
ered for this purpose. The result is the transform 




To - 


m 

P 


To ( qP )Ko ( q) - Ko(qp)lo(q) 


Tq(cl) 


PoO^o* (qPo) " ^KoClPo) " Ko(q)| PoQTo (iPo) “ rolo(qPo) 


(5) 

Inversion of this transform requires knowledge of the zeros of the denom- 
inator in the complex p plane. The pole at p = 0 will contribute the 
steady- state result for the present problem, which is seen to be nonzero by 
physical reasoning. It is 


T p (p,co)-T, 

— - — - = 11m m 
To" T i p-*o 


Io(qp)K 0 (q)-K 0 (qp) I o((l) 


Tq(ol) 


PolKo ' ( iPo ) -mKo ( qPo ) 


-Ko(q) 


PoiTo’ (qPo)- mI o(qPo) 


m 


log ( l/p) 


1 + m log ( l/p 0 ) 


(5a) 


The remaining zeros of the denominator in equation ( 5 ) are real and simple. 
If they are denoted 


q = e +l£ P n > P = "M-n 2 
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and the complex inversion formula 


T c - T i i r 7+i ” P t 
T ° - T 1 ' 2X1 4-i«. 



7 > 0 


is used, then after some algebra, and use of relations for the Bessel functions 
from reference 17 , it is found that 


T c - T 1 = m log(l/p) 

T 0 - T 1 1 + m log(l/p 0 ) 


00 

+ mit 

n=i 


-Mn 2T 


)(^n)|^Po^nJi( Po^n) + m Jo(PoPn) Jo(Phi)Yo(Pri) “ Yo(PP-n) Jo(Pn)j 

p * 2 - ~ 

PoPn J i(Pohi) + xnJ 0 Cp 0 M-n) " U 2 + Po 2 Mn 2 ) J 0 2 (Pn) 


The eigenvalues p n are determined from the equation 

Jp( P~n^ „ JnJ~o(pQpn^ + ^p^n^i^o^n) 

^o(Pn) m -Y 0 (p 0 Pn) Poha^i( Po^n) 


(6) 

(6a) 


From the result for the mixed boundary conditions, case C, given in equa- 
tion (6), the results for cases A and B (see eqs . (4a) and (4b)) can now be 
obtained- First, to obtain case A, m is allowed to become very large, giving, 
in the limit 


t a ■ 


. l0g 4 


T„ - T. 


log (^) 


+ it ; e 
n=i 


" a n • Jo( a n^o(^o a n^ t Jq( P a n)"^o( a n) " ^-n(P a n )Jr>( a n ) J 


'n / o v 


n' o' n' 


Jo (Po a n) “ Jo ( a n) 


(7) 


where 


J o( a n) = J 0 (Po a n) 
^o^n) ^o^o a n) 


and is written for the typical eigenvalue in this case to distinguish it 
from that for case C. 

To determine the solution for case B, write m = h(b/k) and allow the 
product h(T 0 - T 1 ) = m(T 0 - T 1 )/(b/k) to approach unity (corresponding to 
(aF 0 ) =1 in eq. (4b)) while h 0. Then, 
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- = log i + JC 


T - T 

B 

b/k 


00 

l- 

n=i 


-3 n 2 T J 0 (P n )J 1 (p 0 P n )[J 0 (pP n )Y 0 (p n ) - Y 0 ( P p n )J 0 (p n )] 


Po-nC J i (PoPn.) ” Jo (Pn)-l 


(8) 


Jo(Pn) _ J i( PpPn) 

YoOn) ~ Y 1 ( P 0 P n ) 


(8a) 


Determination of Displacement and Stresses 

When the temperature distribution in the material is known, it is possible 
to calculate the displacement u and the stresses a r , cr^- Because of the 
radial symmetry, the shearing stress t t q vanishes everywhere* On this 
account, the maximum shear stress at any point is simply (l # 2)|a r - ag| . 

Plane stress *- If a body is thin, as a plate, in one dimension (say z), 
and all the forces and stresses lie in the plane of the body, a two-dimensional 
problem results, in that all stresses a z , ay z , (? xz vanish and the remainder 
are free of z** (See ref. 1-.) The formulas for displacement and stresses are 
then (ref* l' /l ) 





T 1 )x dx - Ea(T 


T x ) 


E 


1 - v2L 


C 1 (l + v) + C s 


1 - v 


(9c) 


The constants and C 2 in equations ( 9 ) are to be determined by means 

of stress boundary conditions. As mentioned above, these will be such that 
there is a pressure p Q acting at the inner surface r = b, while the outer 
surface, r = a, is to be either stress free or rigidly restrained. For the 
first case, that of free outer surface, the constants are to be determined by 
means of 


11 

1 

V 

O 

at 

r = b 

r r = 0 

at 

r = a 
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Now write 


u 


1 + V 


a(T 0 - T x )a 


i(p^p 0 ;t) + D X P + ~ 


(lOa) 


Or- 


a(T 0 - T 1 )E 


= - — i(p,p q ;t) + — — 


1 - V‘ 


(1 + v)D - (1 - V) 

p 2 . 


(10b) 




f i - T 

7- = •— i(p»p 0 ;t) - ~r + - 1 ■ ■ 

t(To - 'i'jE P 2 l c - l x 1 - v2 


(l + v)D + (l - v) 5| 

P J 


where 


/ » f P I'Cp^t) - T 

I(p.Po;t) = / 1 1 P 3 dp x 

J Po T o - T 1 


(10c) 

(lOd) 


The above boundary conditions give 

D 1 = (1 - v) 


pp^Po* + i(1jP 0 ;t) 
1 - Pn 2 


D 2 = (l + v)p 0 2 


Po* + i ( 1 ^p 0 ;t) 
1 - Po 2 


where 


P 


* = 


o 


Po 

a(T 0 - '1\)E 


If the quantities on the left in equations (lOa) through (lOc) are denoted 
by u*, cr r *, Oq* for the present case of stress-free outer boundary, there 
re suit s 


u* 


1 + v 
P 


l( P } Pq-5 1" ) 


+ 


i(i? p 0 ; t ) 

1 - Po 2 - 


(l - v)p 


+ 


(l + v) 



+ 


Pq 2 Po* ~ 
1 - Po 2 - 


(1 


v)p 


+ 


1 + V 
p 


(Ha) 
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cr r * = 


1 t t \ K1,PoJt) f 

-~2 MPjPoJt) + 2 (1 - 

p2 1 - P r 2 V 


fof ) Pq 2 Pq* 1 . 1_ 
P 2 / 1 - Pq 2 \ P 2 


(lib) 


T - T i . „ i(i,p 0 ;t) t g p 0 ^Pq 


i2n * 


a* = -a* ± + 2 

0 r T o - i - p Q a 


1 - Pr 


(lie) 


Next, the boundary conditions of rigid restraint of the outer surface and 
pressure p Q on the inner surface give 


Di = "(I - v 2 ) Po 


2 Pq* “^(^"■ J Pq-’ 


(l - v) + (l + v)p 


D 2 = (l + v )p r 


(1 - v)p Q * - (1 + v)l(l,p 0 j T ) 


(1 - v) + (l + v)p 2 


Denote now the quantities on the left in equations (lO) by u T , a r 1 , o' . 
Then 


u 


' = ~~ V i (p^Po;t) - (l + v) 


l(l> P 0 ^ T ) 


p w (l - v) + (l + v) 2 


p 2 (1 - v 2 ) Pr 


(l - v) + (l + v)p2 \ p 


I - P 

r\ 1 


(l - v)p + (l + v) -g- 

(12a) 


I _ 


- i i(p,p„;t) + -P 4 v)i(1 ’ p ° ;t) (V 


(l - v) + (l + v)p Q 2 


p * 
lo p o 


(l - v) + (1 + v)p o 2 


(l + v) + 


1 - V 
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Oq' CT r ' - — - 2(1 + v) 

O “ 1 


pQ^'lC U Pp j T ) 


(l - v) + (1 + v)p o 2 


- 2(l + v) 


D 2-n * 

p O F o 


(l - v) + (l + v)p Q 2 


(12c) 


It is interesting to note that when the outer boundary is free, the 
stresses c> r *, o^* do not depend upon the Poisson's ratio v. This is not 
true for cr r ' and Gq 1 , that is, when the outer boundary is rigidly restrained. 
In either case, however, the radial displacement does depend upon the Poisson's 
ratio . 

Plane strain .- Here the idealization is that of a long body in which all 
quantities are independent of axial distance. In this situation, there are 
axial stresses cr z acting at the ends of the cylinder of magnitude 


°z = a r + °0 


These can be removed by a uniform axial stress applied so that the resultant 
force on the ends is zero. By Saint Venant's principle, this self- 
equilibrating system acting on the ends will affect only locally the stresses 
in the cylinder (see ref. 1*0 . 

In order to obtain formulas for plane strain problems from the above 
results for plane stress. It is only necessary to substitute E/(l - v 2 ) , 
v/(l - v), (1 + v)a for E, v, a, respectively. Thus, for example, equa- 
tion (lib) becomes, for plane strain. 


(1 - v) 


I ( P j Pcm ' ) 


I(1,P 0 ;t) 


£of ) , Pq 2 Pq* 

e 2 / 1 - p q 2 



a(T 0 - '1' 1 )E 


The Steady- State Distributions 


In all cases, the steady- state temperature distributions are represented 
by the formula 


T(p,co) - T. 

— ^ 1 = « 1°S i 

i ref 


( 13 ) 


10 



I 


For the three 
Case A 
Case B 
Case C 


heating modes the constants take the following values: 




T ref - T o ” T i i 
T ref = b / k i 
T ref = ■'•o " T i i 


N = 


log 



N = 1 



1 


+ m log 



These temperatures therefore decrease monotonically between p = p 0 and p = 1. 
The stresses, as given in equations (lO), depend upon the integral 


I(P,P, 


>J°°) = / 


p Tfp^oo) - T x 


T. 


ref 


*>l - I 


log P - P 0 2 log p 0 - i(p 2 - P 0 2 ) 


(l4) 

It is interesting to determine whether the steady- state stresses have an extre- 
mum in the interval p < p < 1. 

Consider first the stresses as given in equations (lO) before any stress 
boundary conditions are inq>osed. Then, 


d a. 


aT ref E dp 


js I(p,p o ; “) - ! ** p ♦ fi 


d ° e - -^i(p,p o; ») D = 


aT ref E dp p 


P p 3 1 + v 


An extremum of a r occurs at p = p cr where 


cr 


= P 0 2 (l + log 


+ Tt 


4 Dg 

N i + v 


(15) 


When the value of D 2 found above for a stress-free boundary is used. 




(15a) 
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It is interesting to note that when the internal pressure is neglected, the 
critical point for a r in the free outer boundary problem is completely inde- 
pendent of all elastic and thermal parameters, depending only on the geometri- 
cal quantity p Q , the ratio of inner to outer radii of the body. Also, since 
both p Q * and N are positive quantities (when T re f > 0), it is seen from 
equation (l5a) that a sufficiently high internal pressure p Q can cause the 
position of the extremum, p = p cr , to move beyond p - 1. 


Examination of the derivative of circumferential stress indicates that no 
p in p Q < p < 1 gives an extremum. Hence, the (algebraically) largest and 
smallest values of circumferential stress must occur at the boundaries, at 
least in the steady state . 


In the event that the outer boundary is rigidly restrained, equation (15) 
gives for the critical radius 

2(1 - v)p * 


Per — ^P r 


1 + (1 - v) log (i 


N 


(1 - v) + (l + v)p 2 


(16) 



Figure 1.- Locations of extreme values of 
radial stress in the steady state 
(eqs. ( 15 ) and (l 6 )). 


Again the circumferential steady- state 
stress is monotonic in p Q < P < 1. 
With the rigidity condition, the value 
of Per; the position of an interior 
extremum of radial stress, does depend 
upon Poisson's ratio v, even when the 
internal pressure p Q is ignored. 

In figure 1 P cr is plotted for 
free and for rigid outer boundaries. 
The latter are drawn for three values 
of the Poisson ratio v, and p Q is 
always taken to be zero. It is seen 
that the effect of change of Poisson 
ratio is not great, even over the 
extreme range 0 < v < 0*5- 


Initial Behavior of Solutions 

The series derived above as solutions for the temperatures and their inte- 
grals which are used to determine stresses are not well adapted for use with 
small values of the "time" t, say for t < 1CT 2 . The damping effect of the 
exponential terms is then delayed, and many terms of the series must be evalua- 
ted. Although this is not a serious disadvantage when the summations are cal- 
culated by an electronic computer, it is often useful to have a simple 
expression for the initial behavior of a solution. The applicability of the 
present results is thereby enhanced since cases not covered by the specific 
calcinations given later may be treated. The eventual distributions in the 
steady state are known from the previous section. 
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A particularly advantageous way of obtaining the small r results was 
popularized by Goldstein (ref, 15) . This involves the use of asymptotic expan- 
sions of the Laplace transforms (Goldstein used operational concepts) of the 
various f unctions of interest , followed by term-by-term inversion of the 
results* A discussion of the process, with justification, can be found in ref- 
erence 18. 


Approxima tions for the temperature.- To start, consider the result given 
in equation (5) for the Laplace transform of temperature: 


T C ~ T i _ m 


B 0 (qp) 


where 


T o “ T 1 P mE 0 (qp 0 ) - qP 0 B 1 (qp 0 ) 


Bo(z) = I 0 (z)Ko(q) - K 0 (z)l 0 (q) 
Bj_( z) = I 1 (z)K 0 (q) + K 1 (z)l 0 (q) 


= B 0 T (z) 


( 17 ) 


} (17a) 


The asymptotic formulae for the modified Bessel functions will be needed; from 
reference 17* 


Ic» 


I 1 (w) 

Ko(“) 

K x (oj) 





> 


(18) 


J 


where terms O(l/co 2 ) 
the quantities B Q , B 


in the parentheses have 
x can be found: 


been dropped. 


From these results. 
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B ° ( “ ) - - ^ i (q ' u> t 1 - i (i - V]- e ' u ' a) t 1 + i(i ■ i)] } 
- ife {• (4 '" ) t 1 + 3 (? ♦ i)] + t 1 - § (J + -i)] } 


(l8a) 


J 


The last results can now be substituted in the transform equation (IT) , 
and the following approximate (p -*• oo) f ormula is found: 


T C ~ T x „ _m_ fe-^P-Po) 
T 0 " T x s/PoP L p(q + 7) 


1 1 - p e~^ p ' Po) ~ 

8 p pq(q + 7)- 


( 19 ) 


where 



Po + 8m) 


In this approximation, terms in e -2( l( 1- Po)^ etc., have been dropped, on the 
suppositions that for the present problem there will be most interest in the 
region where p - p Q is small and that 1 - p Q is not small. If the latter 
condition is violated, so that the wall is quite thin, some of the terms 
neglected here could easily be included. These terms represent the influence 
of the outer wall. It is an advantage of this expansion method that such 
influences can be identified and the ones of interest in particular cases exam- 
ined separately. 

The transform in equation (19) can now be inverted; the result is, for 
small values of t. 


T c (p,t) - T x 
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It will "be noted that equation (19) is not fully expanded in powers of 
q "* 1 hut, rather, the denominator (q 4- 7 ) is kept intact. In the limit 
t -+ 0, results of the two processes must, of course, coincide, hut there is a 
physical reason for using the expansion as in equation (19) rather than a com- 
plete expansion in powers of q _1 . In the approximations used above, where the 
effect of the outer wall is neglected, and if curvature is also neglected, the 
problem for the hollow cylinder is like the same problem for a semi-infinite 
slab. The one- dimensional slab problem is expressed by the equations 


5 2 t _ St . 
Sx 2 S T 


T(x,0) = T 1 ; 11m T(x,t) = T. 

^ X->oo ± 


St 

Sx 


= m[T 0 

x=o 


T(0,t)] 


The solution for this problem, in transform notation, is 


T - T, ^-qx 

± = m — 

T o “ T i p(q + m) 


m 


e~q x 

pq 


+ 



1 


m 


e~l x 

pq 



•) 


The partial sums of the last series can be interpreted as successive approxima- 
tions to the solution of the problem in question. For, if 


t(°)(x,t) = T 1 

and a sequence of solutions of the diffusion equation is defined by the condi- 
tions 

T^^(x,0) = T x 
lim T^)(x,t) = T, 

X-^co 1 

- = m[T 0 - T (1 _i) (0,t)] 
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then the Laplace transform of the quantity (T n (x,t) - T 1 )/(T q - T^) coin- 
cides with the nth partial sum of the series expansion above. Thus, there is 
introduced by the expansion in powers of q" 1 an additional approximation in 
which the boundary condition is never exactly satisfied. This would seem to 
explain why the results obtained in the original problem by keeping the term 
(q + 7) _1 are somewhat better than those found by straightforward expansion in 
powers of q -1 , as will appear below when comparison with machine - c omputed 
results is made. For reference purposes, the expansion in terms of q -1 is 


Tp( P, t ) - Tj. 


T, 


T, 


m 

n/PqP 


f 2-r 1|,,2 ierfc 


P - Pq 

2\f~ 




1 

P 


4 1 + gm ) 
Po / 


Ti 2 erfc — — 

2 sfr _ 

(20a) 


The functions i n erfc x are defined as 


i n erfc x 


. n- 

l 


erf c 


Jx 


i dg 


i°erfc x 


erfc 



dg 


Tabulations of this integrated error function can be found in references 1 6 and 
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It is found that the results given by the approximate formula (20) are 
adequate for the smaller values of Biot number m, but when m > 1, accuracy 
falls off. This situation can be remedied by writing 


T 


C " 


T 


o ” 





(20b) 


and expanding the two terms on the right side separately . The expansion for 
T A will be given shortly, and the difference term is approximated by the fol- 
lowing expression: 


T - T 

h _A 

To - T x 


i _ JL (_L + i - i 
87 \ p o 


exp [ 7 ( p - P Q ) + 7 2 t ] erfc ( P g ^_ P ° + 7 V t 


+ 


87 



+ 1 



erfc 


P ~ 

2 -Jt 



(20c) 
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The small t results for cases A and B of heat input are 

Tl *J ® { C 1 + 3 TT 6 ) el,pta(p ■ p ° ) + d2Tlerfc (ttt 2 + a ^) 


. 1 1 ^-P erfc 

^ 8p 2 \Tt 


( 21 ) 


where 


a.^li 

8po 




1 + 1_ 1 - P ) Lrfc p ~ p ° 

b/k Pn/PqP V V 8 p P / l 2J= 


- exp [ 3 ( P - p n ) + P 2 x]erfc P ° + 3*/^ 

2 n/t 


Hf (*) /s exp 


(•° ~ Pq ) : 
4t 


+ (P - P 0 )' ” Q P erfC 

»P 


2*/r 


where 


( 22 ) 


P = 


3 + P c 
8Po 


Comparisons of the above approximate results with some machine-calculated, 
ones appear in figure 2. The latter are based on the series given in equa- 
tions (6), (7) , and (8). The cases m = 0, 0.1., 1.0, 10.0, 00 are shown with 
Po = 0.15. 

It should be mentioned that the quantity plotted as V c differs according 
as m < 1 or > 1. For m > 1, 


and as m -*• 00 


v c - 



T 

T - 


l 


i 


lim V P = V A 

m-*oo 0 A 

which can he verified by conparison of equations (6) and (7). On the other 
hand, for small m. 
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11m 

m->o 


as discussed at equation (8) . 


T c - T l = 
m(T 0 - T x ) 

The quantity 


b/k 

Vq is thus defined as 



when m < 1, so that there is a smooth approach to the limit as m 0. Note 
that the temperature ratio (Tq - T 1 )/(T q - T 1 ) for m < 1 is obtained by 
multiplying the shown quantity Vq by m. This does not hold in the limit 
m = 0, where the quantity graphed is 




_ t b 


T„ 


b/k 


since the tenperature analogous to the ambient T Q is then undefined. Curves 
are drawn for values of p generally near the inner boundary, which is the 
region of greatest interest at small times. The agreement between the approx- 
imate and exact results is seen to be quite good up to about t = 0.01, but the 
approximate results become more and more in error thereafter. 

Also shown, on the portion of figure 2 corresponding to m = 1.0, are the 
asymptotic results from an expansion in powers of q" 1 (eq. (20a)). These are 
seen to be much worse than the ones actually used, which prompted the remarks 
above on choice of expansion for the temperature . 

Although not shown, a few calculations for a disk with a thinner wall 
(p Q = 0-5) indicate that agreement between exact and approximate results for 
small t is even better than for p Q = 0.15. A concomitant which must be 
accepted is that these results will become less adequate for thicker walls, 
that is, for p Q < 0.15 say. If the error becomes unacceptable for such small 
p Q , it may prove profitable to use Goldstein 1 s device (ref. 15 ) in which func- 
tions of (qP Q ) are expanded for small rather than for large argument. 

Approximation for the integral I(p,p 0 ;t).- The integral defined in equa- 
tion (lOd) must be calculated before the displacement and stresses can be 
found. Using the form ( 17 ) for the temperature, it is found that 


[“BoCaPo) - q.Po B i( < lPo) ] X J P i B 0 ( Q.P 3.) ^-P 1 

Po 

_ nj_ pB^qp) - P 0 B 1 (qp 0 ) 
m Bo( < lPo) " qPoB^(qPo) 


(23) 
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Straightforward expansion using the asymptotic formulae (l8a) gives 


i c (p,p 0 ;p) 


m 


pq(q + 7) 


i Dl -q(p-Po) 

sJPo 


m 


8p 2 (q + 7 ) 


1 - ■ 11 + tjM e 


-q(p-po) 


where again 


7 = 7r^-(3 + Pq + 8m) 
8pn 


(24a) 


This result can be inverted to yield 
-^-q( P > Pcv t ) 

1 


= m 


y. 


-1 + 2/ (-^y + e^ 2r erfc (y K fr) 


1 C 1 + k 


87 


+ -3L 1 1 + -p- 


Po. 


- m 


Po l 7‘ 


1 - 2- (l + 1 

87 


-erfc £^Po + 2 y^ ierfc 

2 */r 2^ 


+ exp[r(p - P ) + 7 2 t] erfc 


P - P 
2\[t 


87 


0 + )r L + - 5 =- (l +4) Wi 2 erfc P P ° 


2^t 


(24b) 


The result in equation (24b) is adequate for the smaller values of the 
parameter m shown in figure tut the accuracy is falling off as m 
increases* It is possible to improve the approximation somewhat for larger 
values of m by the same device that was used above for the temperature: 

i c = i A + d c - W (25) 

and expanding the difference expression separately* The expansion of 1^ for 
small t will be given later. The difference expression is, for large q, 

1 qP 0 B i( < lPo) pB x (qp) - PoB x (qP 0 ) 


Ip - 1^ - 


Pq B 0 (qp 0 ) mB 0 (qp 0 ) - qp 0 B 1 (qP 0 ) 


p(q + 7) 


fp. 

I Po 


i+r_(i + 3 + _4 
8q \ P P a 


-q(p-Po) 


- 1 


" 8q (p ? o + 0} 


(25a) 
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The last expression can he inverted to give,, for small t: 

I c - J A ■ p o^ | erfc 1 ^° - exp[r(p - Pq) + 7 £ t ] erfc ^_ P ° + r-s/t) j 

+ 8^J^( 1+ i + 'fe){ exp[r(p ' P o) + 7 2 t] erfc (^r^ 2 + >/t 
- erfc — + 2>s/T ierfc — — -Tsl - i 1 - e^ 2 " 1 " erfc (j-Jr) 

2n/t 2Vt J 7 . 


-i- (-£- + 1 

87 s \Po 


e y2r erfc (t/t) - 1 + 27 (-By 


\l/2 


(2?b) 


In order to use this expression 
to have the small t expansion for 
that for Ig. 



in the calculation of Iq, it is necessary 
1^. This is given below, together with 



2t ly/ 2 ierf c 


P - Pp 

2n/t 


+ 


1 

2 



+ iL \ -ri 2 erfc — ^2 

P/ 2^ J J 


( 26 ) 


i b(p>Po> t ) 



4ri 2 erfc — — 

2\fr 


3 




t 3/2 i 3 erf c 

2-Tt 


(27) 


Finally, it is necessary to know the values taken by the integral 
I(p,P q ;t) when p = 1. That it is not strictly correct simply to put p = 1 
in the foregoing approximate expressions is obvious because the terms in 
-q ( i-p) 

e were eliminated on the basis that p is not close to unity. Thus, 

by the Wronskian of I ( z ) , K Q ( z ) 

Bi(l) = q." 1 


while the approximation (l8a) gives, when the negative exponential is 
neglected , 


B x (q) 
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Thus, by putting p = 1 before using the asymptotic expressions, it is found 
that 

1/2 


i a (i, Po ;t) - Po + 2 p 0 " 

+ 2 - f) Ti 2 erfc — — —^2 \ 

\Po J 2*>/t J J 


2 t 1/2 erfc 1 ~ P Q 
2 «/t 


(28a) 


I b (i,P q ;t) = t - -|= 


^ L 


4 ri 2 erfc — — 

2 -Jt 


-2- + 1 r ' i erfc 


3 /2,3„ „ 1 - Po 


2n/t 


(28b) 


I 0 (i,p o ;t) <.{i 


1 ‘ 87 \Po + 1 


e^erfc (tt i/2 ^ -1 + 27 


+ ± (JL + i) T l- 

87 \Po / J 21/2 
r Pn 


7(i-P c )+7 2 t ^ fl - Po 1/2 

e erfc ( — + 7 t 

2t 1/2 


erfc — — + 27T 1 / 2 ierf c — — 

2 t i/2 2t 1/2 


(28c) 


I C _ X A| n _ 


Po 


P=i 7 


1 - A. (-L + 1 

87 \Po 


1 - e^erfc 7 '•/?) - ^ (i + l) (| 


l/2 


7 \ 27 p 0 . 


' 7(i-p 0 )+7 2 t /l - po 

e erfc ( ° 


+ 7*/t > ) - erfc - £ 

2n/t / 2^ J 


Vt ierfc 1 ~ P ° 

7*/Pq 2 n/t 


(28d) 


Figure 3 shows comparison of exact and approximate results for the inte- 
gral I(p,p 0 ;t). They appear to be of better reliability than those shown 
earlier for the tenperature , probably because they represent integrated quan- 
tities. The same convention on multiplication by Biot number m applies here 
as for the temperatures in figure 2. 


CALCULATION OF TEMPERATURE AND STRESSES 


The calculations amount essentially to summing the series derived in pre- 
vious sections, such as equation (6), and determining the integral I(p,p q ;t) 
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(see eq. (lOd)). A series expansion of this integral will he given below. 
Although the series are convergent, the favorable damping due to the factor 

e in each term is somewhat delayed when t is small, say t < 10” 2 . In 

this event it is sometimes necessary to consider upwards of 30 terms of the 
series in order to achieve four significant figures- Fortunately, the terms 
alternate in sign so the maximum error (and its sign) caused by stopping at a 
given term is known. 


Calculation of Eigenvalues 

Determination of the sums indicated in equations (6), (7)> (8) and of the 
integral in (lOd) requires that the eigenvalues M- n ( a n^ P n ) found. These 

are defined by equation (6a) ((7a), (8a)). It is possible to derive approxi- 
mate expressions for the eigenvalues when n is large by making use of the 
asymptotic expressions for the Bessel functions (see ref. 17) • Denote, for 
large values of the argument, 


y/2 

J 0 (x) ~ [•—) [P Q (x)cos Xj. - Q 0 (x) sin xj 
\l/2 


jfx) ~ VjxxJ [P 1 (x)sin x x + Q 1 (x)cos xj 


\ l/ 2 


Y 0 (x) - (jLJ [P Q (x)sin x ± + Q q (x)cos xj 
Y-Jx) ~ (lx) [-P 1 U)cos x 2 + Q i (x)sin x i ] 






(29) 


where 


P (x) = V 

° ^ (2X) 2111 


m=o 


Q • (x) =V (-J-) m (j,2m + 1) 


m=o 


2 2k kl 


(j>k) = — i " l2 )(^j 2 - 3 2 ) • • • [kj 2 - (2k - l) 2 ] 


and 



2k 



These expressions are now to be inserted in equation ( 6 a). The result is 


tan(l-p Q )x 


m[P 0 (x)Qo(Pox) - P 0 (p 0 x)Qo(x) ] - P 0 x[P 0 (x)P ;i (pox) + Qq(x)Q 1 ( p p x) ] 
m[P 0 (x)P 0 (p 0 x) + Q 0 (x)Q 0 (p 0 x) ] + p 0 x[P 0 (x)Q ;L (p 0 x) - P 1 (p 0 x)Q 0 (x) ] 


( 30 ) 


where x now stands for the eigenvalues \± n > To solve for those eigenvalues, 
a process of successive approximation can be used, starting by neglecting all 
but lowest powers of x -1 in the expression on the right# For example, in the 
first approximation 


tan(l - p Q )x ~ - j 

This relation shows the characteristic behavior of the well-known transcenden- 
tal equation 


tan x = -x 

That is, the nth root approaches [n - (l/ 2 )]jr as n increases. In the 
present case, the first approximation to the eigenvalue M-n for large n is 



•with n, and with Biot number m 
(eq. (6a) , p Q = 0 . 15 ) . 



Further approximations can be obtained 
if desired, if the result (31) is sub- 
stituted into equation (30) and terms 
involving higher powers of 

|[n - (1/2) | are determined. 

Such approximations were not used 
(except for checking purposes) in the 
present calculations. It was found as 
easy to use the complete frequency 
equation ( 6 a) and determine the roots 
by means of a simple interpolatory pro- 
cedure . A standard program for Bessel 
functions was used when the argument 
x < 20 , and extensions made by use of 
the asynptotic formulae (29) * Figure k 
shows the variation of the eigenvalues 

(actually M* n /[n - (1/2)^]) with n. 
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and also with the Biot number m which appears as a parameter in the fre- 
quency equation. Limiting values m -> co where p n -+ Oq and m = 0 where 
M- n = P n are The eigenvalues for m = 100, 0.01 are nearly indistin- 

guishable from those for m -> oo, 0 , respectively# This particular graph is 
drawn for p Q = 0.15, so the limiting value as n oo is 


lim 

n-K)o 




1-18 . . . 


Summation of Series Expressions 

Once the eigenvalues are known,, one can continue with the summation pro- 
cess. In this, it was found adequate, for t > 5 X 10 ~ 4 .? to provide for using, 
at most, 50 terms of the series. The computation halted automatically at a 
value of n such that the corresponding term was less in absolute magnitude 
than 1 X 10 " 7 - 

The integral I(p,p 0 ;t) defined in equation (lOd) can also be expressed 
as an infinite series. Using equation (6) for temperature, it is found that 

IqCp^PoJt) 





00 O 

v ”^n T 

- 2m 2 \ -§ 

n=i ^ tmJ 0 (p 0 ^ n ) + P 0 ^ n J 1 (p 0 ^ n )] 2 - (m 2 + P 0 2 ^ n 2 )J 0 2 (^ n ) 


+ Jtmp 


00 

1 

n=i 


J 0 (n n )[mJ 0 (p 0 ^ n ) +P 0 ti n J 1 (p 0 P rl )][J.,(p^ n )Y 0 (^ n )- Y 1 (p l i n )J Q (p n )] 
[mJ 0 (p 0 P n ) + Po^Po^rh 2 ' + 


( 32 ) 

The results for cases A and B can he obtained from equation (32) by suitable 
limiting processes performed on m. Results of the calculations outlined above 
will be presented in the next section. 


RESULTS OR COMPUTATIONS 


There are an enormous number of possibilities for calculating the 
stresses, even though the parameters E, a(AT), a, k appear in dimensionless 
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ratios. A selection from the remaining ones* m, v, p Q , was made which was 
thought to give a representative selection of cases for different regimes. It 
is hoped that the cases detailed here will enable one to obtain at least a 
qualitative notion of the magnitudes and distributions of temperatures, 
stresses^ or displacements in a wide variety of situations. There are available 
also the fairly simple results for the steady- state values and for the initial 
variations derived above. These formulae are of course quite general insofar 
as assigning the parameters m, v, p Q is concerned. 


Effect of Variation of Biot Number m 


Free outer boundary.- The first case examined is for a fixed geometry 

(p Q = O. 15 ) and Poisson 1 s ratio (v = 0.25)* with the Biot number m allowed 
to vary between its extremes of 0 and 00 . All quantities of interest are shown 
in figures 5 through 9 , including curves of the integral I(p,p q ;t) (see 
eqs. (lOd) and (32)). The curves are drawn so that the variation of each quan- 
tity with p is shown for values of t covering the range 0.0005 < t < 2. 

The latter value of t is one for which the system has essentially reached its 
steady state, as indicated by the near identity of the calculations for t = 1 
and t = 2. The integral I(p,P q ;t) is shown so that the present calculations 
can easily be extended to stress boundary conditions other than the vanishing 
of radial stress at p = p Q and p = 1. (Some results for the boundary condi- 
tion of a rigid outer wall will be shown later.) Also, the internal pressure 
p Q is zero in the calculations for the stresses and displacement; its effect 
may be included by use of equations (ll) or (l2) . 


The values of m for which results are shown are m = 0(.0l), 0.1, 1, 10, 
oo(l00) . The numbers in parentheses indicate that the corresponding results 
differ so little from the results for m = 0 and m -► 00 that they are plotted 
on the same graph, insofar as they are distinguishable. As discussed above, a 
difference exists in the dimensionless groups representing temperature, dis- 
placement, and stress according as m < 1 or m > 1. For m < 1, m is included 
in the denominators in order that the dimensionless quantities V c approach 
the proper values at m 0, corresponding to the constant flux case. For 
m > 1, m is not included in the dimensionless groups so that the proper limits 
result as m 00 , corresponding to the case of sudden change of temperature at 
p = p Q . Thus, for example, 

v c = m " 1 ? s — Jr for “i 1 

o 



for m > 1 


This also applies to l(p,p 0 jr). 

The graphs of the radial temperature distributions show characteristic 
behavior; the temperature at first rises quickly in a thin layer of material at 
the inner surface, and the tengperature gradient across this layer increases 
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■with m. The curves of displacement show, for all values of m, a rather weak 
maximum at an interior point which moves toward the outer boundary as time 
increases. The radial stress, which is here zero at both boundaries because 
p Q = 0, shows similar behavior, although the extremum is much accentuated as 
t increases. It is, for T Q > T-j_, compressive throughout, and its position of 
numerical maximum moves from the inner wall for t « 1 toward the value 
p = p cr predicted in the section on steady-state behavior (eq. (l5a)). 

In the case of the circumferential stress, a flat extremum occurs at first 
near the inner boundary, moving to and past the outer boundary as time 
increases. The value at p = p Q is the largest compressive stress, and the 
largest tensile stress occurs eventually at p = 1. The maximum compressive 
stress is larger in magnitude than the largest tensile stress, and also larger 
by a factor of 3 to 5 (depending on m) than the maximum radial conpressive 
stress. The behavior of the circumferential stress distribution is in some 
respects similar to that of the temperature distribution for small values of 
t. That is, it is nonvanishing only in a thin layer near the inner surface, 
the magnitude at the inner surface and the thickness of the layer increasing 
with time. For the more intense heating (m > l) the maximum conpressive cir- 
cumferential stress, at p = p Q , no longer occurs in the steady state as it 
does for the m < 1 cases. This behavior carries over into the stress differ- 
ence (cq,* - cr fl *) , and figure 10 shows this maximum value as a function of time 

for the values m = 1, 10, oo- The 
curve for m = 1 is characteristic of 
those for the smaller values of m 
also. 

Finally, in figures 5 through 9 j> 
the stress difference (a r * - cr^*) 
gives a measure of the maximum shearing 
stress to be found at any point. At 
the boundaries, this quantity is domi- 
nated by the circumferential stress, 
and it falls to zero and changes sign 
somewhat closer to the inner boundary 
than does 

Rigidly restrained outer boundary.- 

In order to gain some notion of the 
effect of a rigid outer boundary, a few 
curves will be given in which the 
results for a free boundary and a fixed 
boundary are compared. These are drawn 
only for m = 0, 1, oo and for values 
of t which show the starting behavior and steady state as well as the transi- 
tional interval. Also, only the displacement and stress difference are shown. 
Figure 11, where the displacements for the three values of m are given, shows 
that initially (t = 0 . 005 ), the differences between the displacements in the 
cases of zero and conplete restraint are rather small. However, as time goes 
on, the differences increase and, as expected, the free displacements are ulti- 
mately considerably greater throughout than those for a rigid restraint. These 



Figure 10.- Behavior of maximum shearing 

stress (a r * - a 0 *)p_ p tor several values 

of m(p Q = 0.15, v s O. 25 ). 


33 


.1 



Figure 11.- Comparison of displacements u* 
and u* for free and fixed outer boundary , 
respectively (p 0 = 0.15, v = O.25). 



Figure 12.- Comparison of maximum shearing 
stresses (a r * - cr^*) and (o J . 1 - Gq' ) fo: 
free and fixed outer boundary, respec- 
tively (p = 0.15, v = O.25). 
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effects increase in magnitude with the intensity of heating, as would also he 
expected. 

Figure 12 illustrates the differences between the values of (a r ' - Oq ' ) 
and (cr r * - cT0*), corresponding to fixed and free outer boundaries, respectively. 
The variation is seen to be similar over the range p Q < p < 1, the fixed 
boundary giving the higher magnitudes. Again, differences are heightened by 
increasing t and by intensification of heating, or increasing the Biot 
number m. 


Effect of Poisson's Ratio v 

In the free-boundary case, only the displacement u* is affected by a 
change of Poisson's ratio. Figure 13 shows the magnitude of the effect for the 
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Figure 13*“ Effect of Poisson 1 s ratio v 



extreme cases m = 0, in all cases, the magnitude of the displacement 

increases with v. In the most severe heating (m = 00 ) , the change in displace- 
ment due to a change in v is sizable at its maximum, even for the smallest 
value of t (t = 0.005) considered. 

Similar results for the stress difference (cr r * - cr^*) in the fixed bound- 
ary situation appear in figure l4. Here, for t = 0.005 and m = 00, the effect 
of altering v does not appear to be proportionately as great as it was for 
the displacement results shown above. 


The Influence of Wall Thickness 

The differences between the solutions under various heating inputs and 
restraint conditions for differing values of the geometrical parameter p Q are 
somewhat more difficult to predict than those found in the parameter surveys 
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considered above. First, figures 15 and l 6 show stress differences for widely 
different wall thicknesses (p Q = 0.05 and O.90) for a free outer boundary, and 
for the extreme heating conditions m = 0, 00. In the case of the thin wall 
(p Q = 0.90, fig. 15); the steady- state distributions appear nearly linear for 
both heating conditions, while the starting distributions (t = 0.0005) are 
somewhat simpler than those given previously for p Q = 0.15; that is, no trace 
of an internal extremum in stress difference remains. 

In figure 1 6, the results for stress difference in a very thick ring are 
given. Here the stress gradient near the inner surface is much greater than in 
the thin-ring of figure 15 - A definite interior extremum exists for small t, 
and there is a greater concentration of the stress near the inner surface at 
the initiation of heating. 

Figures 17 and l8 show the quantities corresponding to those in figures 15 
and 16, only for the stress-boundary condition of rigid restraint at p = 1. 

The character of these figures is analogous to -that of the corresponding ones 
for a free outer boundary, with the expected larger values at p = 1 due to the 
rigid restraint. 
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CONCLUDING REMARKS 


The present calculations should be useful in determining heating which 
vi 11 keep the stresses below some allowable maximum,, as illustrated for other 
configurations by Heisler in reference 1. 

The thermal stresses can become considerable , as the following example 
will show. For the case shown in figure 7* the maximum stress is given by the 
value of oq * at the inner wall in the steady state. It is 

|<W| = (0.5)a(T o - T X )E 


For several materials, such as aluminum or steel, the expansion coefficient a 
has a value roughly 1XL0~ 5 /°K to 2X10" 5 /°K. If the temperature difference is 
T Q - T 1 = 200° K, the maximum stress is in the range EX10" 3 to 2EX10" 3 . (if 
the object considered is a pipe rather than a disk, this stress is divided by 
1 - v. For the materials mentioned, this amounts to taking about 4/3 the 
ab o ve re suit s • ) 
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Now, as discussed "by Nadai in reference 13 (ch. 13), the various factors of 
dependence of material constants upon temperature and the onset of plastic flow 
should be taken into account when thermal stresses become severe. If it is 
desired to keep the stresses below the level at which these new effects enter , 
the calculations above will serve to define limits on the heating, or dimen- 
sions, or other parameters which will so confine these thermal stresses. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Sept. 17, 1965 
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APPENDIX 


B o (z) 

B i(z) 
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J o > Y o 
k 
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Po> P 0 * 
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t 


t a> t b' t c 

T o 

T i 

u 


V, 
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SYMBOLS 

outer and inner radii of annular disk or hollow cylinder 
I 0 (z)K 0 (q) - I 0 (q)K 0 (z) 

I i (z)K 0 (q) + I^qjK^z) 

(p 0 - l)/p G (eq. (21)) 

Young ' s modulus 

surface heat-transfer coefficient (eq. (4c)) 

integral arising in stress calculation (eq. (lOd)) 

modified Bessel functions of zero order, of first and second 
kinds 

Bessel functions of zero order of first and second kinds 

thermal conductivity 

Biot number, hb/k 

Laplace transform variable 

internal pressure and dimensionless internal pressure 

pl/2 

radial coordinate 
time 

temperature 

boundary conditions, equations (4) 

temperature of medium in r < b 

initial temperature of solid 

displacement in thermally strained body 

dimensionless temperature ratio, (T^ - T i )/(T Q - T ) 


4o 


V B 

Vc 

Z 

a 


a n 

3 

Pn 

7 

K 

v 


P 


Po 

°r > °Q 


( )* 
( )' 


(t b - T^)/ (b/k) 

(T c - T 1 )/(T 0 - T x ) if m > 1; (T c - T 1 )/m(T 0 - T 0 ) if 

axial distance along a cylinder 

coefficient of thermal expansion 

eigenvalue, equation (7a) 

(3 + p 0 )/p 0 (eq. (22)) 

eigenvalue , equation (8a) 

o^-(3 + P 0 + 8m) (see eq. (19) ) 

8p 0 

thermal diffusivity 
eigenvalue (eq. (6a)) 

Poisson's ratio 

dimensionless radial coordinate ; r /a 
b/a 

radial and circumferential stresses 
axial stress 

dimensionless time variable; rct/a 2 

stress condition of free outer boundary 

stress condition of rigidly restrained outer boundary 


m < 1 


4l 
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